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Abstract 

A new approach for the design and control of flexible space structures is described. 
The approach integrates the structure and controller design processes thereby providing 
extra opportunities for avoiding some of the disastrous effects of control-structures inter- 
action and for discovering new, unexpected avenues for future structural design. A linear 
control formulation based on Boyd’s implementation of Youla parameterization is 
employed. Control design parameters are coupled with structural design variables to pro- 
duce a set of integrated design variables whose values are selected through optimization- 
based methodology. A performance index reflecting spacecraft mission goals is formu- 
lated and optimized with respect to the integrated design variables. Initial studies have 
been concerned with achieving mission requirements with a lighter, more flexible space 
structure. Details of the formulation of the integrated design approach are presented and 
results are given from a study involving the integrated redesign of a flexible geostationary 
platform. 
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Introduction 

The “Integrated Control/Structure Design” study was performed as Task #3 of the 
NASA contract “Aircraft and Spacecraft Guidance and Control Technology”, No. NAS1- 
18762. The objective of this effort was the implementation and demonstration of a pilot 
computer code for the integrated design of actively controlled structures. Recognized 
from the outset was the need for a software environment that would facilitate the multi- 
disciplinary analysis involved. The Integrated Analysis Capability (IAC) software [1] 
was chosen to provide this framework. The IAC utilities provide not only an efficient 
data base system for engineering analysis but also easy access to analysis modules such 
as model generation, finite element analysis, control design, and simulation tools. 

To accomplish the stated objectives, the following tasks were entailed (as written in 
the statement of work): 

1. selection of a representative pilot integrated design problem and the appropriate 
analytical methods and computer software, 

2. development and enhancement of the necessary IAC interfaces for the pilot inte- 
grated design problem, 

3. demonstration of the pilot computer code for the selected design problem, 

4. study of the issues involved in the formulation of practical performance functions 
and constraints, and 

5. documentation of the work performed and plans for follow-on efforts. 

Several significant accomplishments were made in the course of this research. By 
developing a method for calculating the open-loop modal response sensitivities without 
differentiating the structure’s eigenvectors, we were able to avoid the difficult and 
involved method that is normally employed. By using a control formulation (based on 
the Youla parameterization) in which the controller responses themselves are parameters, 
we implemented a promising new approach and verified its usefulness for the kind of 
problems considered here. By implementing a software tool that can not only perform the 
integrated analysis of an actively controlled structure but also optimize its design for a 
wide range of possible objectives and constraints, we demonstrated that such optimization 
is practicable and that significantly better designs can result. 

Besides meeting the original objectives, we have learned much about the problem 
and have identified several areas that may warrant further investigation. This document 
is the final report of the task and describes our objectives, methods, and findings. 
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Analytical Approach and Implementation 

The usual approach for designing large, flexible, space structures with active vibra- 
tion suppression is a two-step procedure. The first step is to determine the structural 
design based on structural costs and open-loop dynamic behavior. Then, given this 
structural design, the second step is to design the control system to obtain the desired 
closed-loop dynamic behavior. These two steps may be iterated until a suitable design is 
obtained. 

This sequential approach is fairly straightforward; however, it makes it difficult (and 
entirely reliant on the designer’s experience and intuition) to exploit any coupling that 
may naturally exist between the two design problems. Because the structural design vari- 
ables and the control design variables are considered separately, information regarding 
this coupling is not available to the designer. Consequently, there is no reason to expect 
this approach to converge even to a locally optimal design for most constrained problems. 

The poor performance of the sequential approach is easily demonstrated. The simple 
optimization problem illustrated in Figure 1 has only two design variables and a single 
constraint; its feasible space lies above the constraint line shown. We might imagine that 
x is a structural design variable and that y is a control design variable. Consider the sit- 
uation shown: the design currently under study lies on the constraint at the point A. Now 
if we are restricted to moving only along the x or y axes, as would be the case if the 
sequential approach were used, then the design cannot be improved without violating the 
constraint. A clear path to the minimum exists, but it can be taken only if the design vari- 
ables can be modified simultaneously. 
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Figure 1: A simple optimization example. 


For this reason, it seems likely that improved designs could often be found by con- 
sidering simultaneously the complete set of design variables. This has been apparent to 
many, and a large number of studies are described in the literature. Some of this work 
has been directed toward formulating the combined problem in a manner similar to those 
of modem optimal control, several examples of which arc described in References [2-5]. 
Unfortunately, these approaches tend to be quite restrictive in terms of the kinds of prob- 
lems to which they apply. These restrictions typically preclude problems of practical 
complexity. 

An alternative approach to simultaneous optimization is to use the methods of 
numerical optimization. This approach has a distinct advantage in that it is much less 
restrictive on the kinds of design variables, objectives, and constraints that can be consid- 
ered. References [6,7] describe several applications of this approach. 
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Attempts to use numerical optimization methods, however, have been hampered by 
several difficulties. Because the objective and constraints are typically defined in terms 
of the closed-loop response of the structure, a complete structural and control analysis 
(with sensitivities) is required for every trial design. The effort involved can be pro- 
hibitive. This is particularly relevant to this multidisciplinary problem, since integrated 
analytical tools are generally unavailable. This is exacerbated when numerical differen- 
tiation is used to obtain the necessary gradients. Moreover, the dimensions of the design 
space and the complexity of the objective surface may also frustrate these optimization 
attempts. 

The approach taken here for the integrated design problem is shown graphically in 
Figure 2. This approach was implemented in a computer program called COSTAR. 
Using a nonlinear programming method, an optimization module performs a search of 
design space, attempting to minimize the specified objective function while satisfying 
various constraints on its design and performance. The design space is comprised of both 
structural and control variables, and the optimizer is free to modify them simultaneously. 
We have attempted to avoid the principle obstacles to this type of approach by using 
efficient, integrated analysis tools, analytical gradient calculations, and design variables 
that tend to reduce the complexity of the objective surface. 
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Figure 2: The COSTAR conceptual layout 

The functions that describe the objective and constraints are defined a priori by the 
user. Although this explicit description of the problem may require a significant effort 
from the user, the hallmark of this approach is that the type and form of the objective and 
constraint functions is restricted only in that they must be continuous functions of the 
design variables. In the demonstration problem, functions such as peak error and peak 
torque” were included along with more conventional performance functions such as “total 
mass” and “mean-square error”. 

Structural Modeling 

In the COSTAR code, equations of motion for the structure are obtained by using the 
finite element method (FEM). A description of the structure is used by the FEM code to 
construct mass and stiffness matrices that characterize the structure. This description 
takes the form of several matrices that contain the connectivity, geometry, element prop- 
erties, and material properties. 
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Table 1 lists the matrices that combine to describe the structure. As indicated, sev- 
eral of the matrices are functions of the design variables vector v. Thus different choices 
for the design variable values result in different structures. Throughout this report, 
derivatives with respect to the design variables are indicated by a prime. 


Matrix 

HB3SSESBI 

Description 

X 

y/_ 3X 

Ai - dvt 

Locations of the nodes 

P 

¥- 
‘ ' dVi 

Element properties 

E 

r /_ M. 

Ei — 

a V,- 

Material properties 

A 

, dA 

AiS dt 

Element orientations 

K 

|£-o 

dvi 

Element connectivity 

P 

dvi 

Element property pointers 

X 

f*-0 

dvt 

Material property pointers 

X 

p-m 0 

dvi 

Element types 


Table 1: Structural Definition Matrices. 


These matrices are themselves generated through the use of PATRAN [8] — a finite 
element pre-processor. The input to PATRAN (its “session file”) is kept in a parametric 
form so that changes to the values of the design variables can be easily accommodated. 
(The MACSYMA symbolic mathematics program [9] is used to update the PATRAN 
input file.) For reasons of efficiency, however, this kind of pre-processor cannot be used 
within the optimization procedure. The overhead involved in running this code would 
add significantly to the total execution time. Moreover, a finite difference procedure 
would be required to obtain sensitivities, adding greatly to the execution time and to the 
numerical errors in the sensitivities. 

In many previous approaches to the simultaneous optimization problem, this problem 
has been avoided by representing the mass and stiffness matrices as explicit functions of 
the design variables. For example, it is common practice to prescribe them to be linear 


8 






functions of the design variables. Although this simplifies the analytical task, it imposes 
severe restrictions on the kinds of design variables that can be considered. In particular, it 
precludes the use of design variables that control the shape of the structure. 

The COSTAR implementation avoids these difficulties by representing the structural 
model (that is, the description matrices listed in Table 1) as a linear combination of a set 
of model variables. The nodal locations, for example, are expressed as 

* = * + + H' 1 ' 12 -'' 2 * + • (1 > 

where /i is the vector of model variables and the overbars refer to the baseline values. By 
specifying the model variables as possibly nonlinear functions of the design variables, 

Mi = Mi(v) ( 2 ) 

relatively few restrictions are placed on the models. With this approach, the modeling 
software is used only to generate the model description. Updates for subsequent analyses 
(such as for optimization) proceed rapidly with only simple matrix operations involved. 
The sensitivities are obtained as 

* + Hw- ! + - (3) 

and similarly for the other definition matrices. 

Finite Element Analysis 

Once the model description is obtained, the next step in the integrated analysis is to 
generate the equations of motion for the structure. These equations of motion are given 
by 

MX + Cx + Kx = F (4) 

The task of the structural analysis module is to assemble the mass, damping, and stiffness 
matrices as well as their sensitivities. While computer codes for finite element analysis 
abound, few are designed to calculate sensitivities and fewer still can calculate the sensi- 
tivities analytically. Because of the importance of fast, accurate sensitivity calculations 
for the integrated design problem, we chose to develop our own special purpose structural 
analysis code (called SSA). 

The SSA routine has only beam and concentrated mass elements, but it can provide 
analytical sensitivities of the mass and stiffness matrices for any combination of continu- 
ous design variables, such as nodal locations, element properties, and material properties. 
A lumped mass formulation is employed (resulting in a diagonal mass matrix), and the 
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stiffness matrix is stored in a banded form. Using the sensitivities of the definition matri- 
ces as calculated in Equation (3), the sensitivity of the stiffness matrix with respect to the 
design variable v, is 


K '. = y y gg 

' ' j 7 ^* 1 * 




dK 


« + 1 S d[P]jk 


[ft ]/* 


(5) 


and the mass matrix sensitivities are obtained in the same way. 

In COSTAR, the damping matrix C is never assembled. Instead, the common modal 
damping approach is used wherein the damping ratio of each mode is specified directly. 
This implies that the eigenvectors of the undamped problem also diagonalize the damped 
problem, which is the case when the damping matrix can be expressed in terms of the 
mass and stiffness as 


C^M^akM-'K) 1 (6) 

i 

with arbitrary scalar coefficients a, and for any integer values of i. Rayleigh damping is a 
special case of Equation (6), with /= { 0, 1 }. With modal damping, the equations of motion 
become uncoupled and their solution is greatly simplified. 

Usually, it is enough to just specify the modal damping ratios and then solve the 
uncoupled system. In COSTAR, however, it is slightly more complicated. As described 
below, the response sensitivities are obtained in modal form but without requiring deriva- 
tives of the eigenvectors. However, this requires the use of C', the sensitivities of the 
damping matrix. Thus Equation (6) must be differentiated and used to calculate C'. In 
the demonstration problem, j=l was used in Equation (6); this results in damping ratios 
that are proportional to frequency. 


Eigenvalue Extraction 

As mentioned above, the open-loop equations of motion are transformed into modal 
coordinates (with mass normalization) as 

qi + 2 Zmqi + cofqi = <p[F (7) 

This transformation both uncouples the problem and enables model reduction through 
modal truncation. At the same time, eigenvalue extraction is often the most computa- 
tionally expensive part of the open-loop analysis task. Because this step must be per- 
formed repeatedly for each new design, the efficiency with which it is performed can 
determine the practicability of searching for an optimal design. For this reason, we have 
attempted to perform the eigenvalue extraction as efficiently as possible. To do this, we 
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have taken advantage of the fact that the eigenvectors from a previous design are usually 
good approximations to those of the current design. 

The existing EIGEN module within IAC is a very efficient eigenvalue extraction 
code. Nevertheless, its Lanczos algorithm does not benefit from good a priori estimates 
of the eigenvectors. Instead, we have implemented a subspace iteration approach (called 
SSI) that makes successive improvements to a set of starting vectors to eventually con- 
verge to the true eigenvectors. The set of equations 

K^k + 1 =M®k 

1 = 

Mk+\ = <Pfc+iM<2W (8) 

'K.k+iWk.+x - Mt+\ y/jr+iAt+i 
i = ^fc+i Yk+i 

are solved for successive values of k until converges to the true matrix of eigenvectors. 
By using the eigenvectors from the previous design as starting vectors for the present 
design, this convergence is often achieved in only a few iterations. 


Open-Loop Responses 

In the COSTAR implementation, the closed-loop transfer functions needed for the 
objective and constraint evaluations are obtained directly from open-loop transfer func- 
tions of the structure via the so-called (2 - dcsign approach described in the next section. 
One advantage of this approach is that it does not require the assembly of a state- space 
(first order) model of the structure. Since the open-loop equations of motion have been 
transformed into modal coordinates, the open-loop transfer functions can be calculated 
very efficiently. 

In addition to the transfer functions, we also need their sensitivities with respect to 
the design variables in order to compute the sensitivities of the objective and constraint 
functions. To determine these sensitivities, it has been common practice to differentiate 
the modal equations of motion (Equation 7) to form an additional set of equations 

q{ + 2CiCO,qi + cofqf = (ftV)' - *)'$; - 2(O t (o'q t (9) 

that characterize the modal response sensitivities Qi . This system has the same frequen- 
cies and damping ratios as the original system, but the right-hand side contains not only 
the forcing function sensitivities but also the responses obtained from solving the original 
system. In this way, the response sensitivities can be obtained from 
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( 10 ) 


x' = O' q + Oq' 

Unfortunately, there are well-known difficulties associated with the calculation of 
eigenvector derivatives [10]. These difficulties are manifested when repeated eigenvalues 
exist — a common condition in practice. 

For COSTAR, we have implemented an approach that appears to offer significant 
advantages. If the physical coordinates of the original problem in Equation (4) are trans- 
formed by some arbitrary constant matrix ©, so that 


x = 07] (11) 

then the equations of motion may be written as 

0 T M0ij + e r C0T] + 0 r K0ri = 0 t F (12) 

By differentiating Equation (12) with respect to the design variables, keeping in mind that 
the matrix © is constant, we obtain 

0 T M07i'+ 0 T C0i]' + 0 t K0ji' = 0 T (F'-M'0ij-C'0T}-K'0n) (13) 

We have not specified the transformation matrix 0; let us now consider it equal to 
the eigenvectors of the original system. With this choice, we can see that Equation (13) 
reduces to 

nV + 2Ci(0im' + cofvi = - m'&v - c'oi] - k'o 77) (14) 

and the sensitivities of the physical responses are 

x' = Ot }' (15) 

The coordinates 77 are equal to the modal coordinates q , but their sensitivities are differ- 

ent. 

Comparing Equations (9) and (10) with Equations (14) and (15), the advantages of 
the latter approach are obvious. In both cases, a second set of equations must be solved in 
which the original responses appear on the right-hand sides. With the second approach, 
however, the physical response sensitivities are obtained without computing eigenvector 
derivatives. Based on our limited experience, this approach seems to work well even 
when O contains a truncated set of mode shapes. 
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Controls Analysis 

A recently developed method for controls synthesis, that is here called Q-design, is 
based on the (^-parameterization (or stable factorization) theory developed by Youla [11]. 
With this parameterization, all possible stabilizing controllers and input/output maps can 
be expressed as functions of a stable parameter Q. With this formulation, the closed-loop 
transfer functions are affine in Q — an important property for our purposes. 

For a given plant P yu and nominal controller K nom , the stable, co-prime factorization 


P yu =ND ~ 1 =D~ l N 

k nom = r- 1 x=xf - 1 

may be chosen such that 


Y 

X 

D 

-X ' 

. -N 

D 

N 

Y 


In that case, the set of all stabilizing controllers is given by 


(16) 


(17) 


K = {(r~ oft) _1 (x +QD)\ Q Stable) (18) 

Furthermore, all achievable closed-loop input/output maps are described in terms of the 
affine Q parameter via 


H Z „ = T X +T 2 QT 3 (19) 

where the nominal system is characterized by 

T i — P zw~P zuPXP yw ~ H 2 w \q=0 

T 2 = -P zu D=H zv (20) 

T 3 = DPyw = Hew 

The Q parameter may be thought of as a tuning parameter for the nominal controller. 
Figure 3 depicts its connection to the nominal controller. Here P represents the plant and 
ATnom represents the nominal controller. The plant has control inputs u, exogenous inputs 
w, regulated outputs z, and measured outputs y. The exogenous inputs include signals 
such as commands, disturbances, and sensor noise. The regulated outputs are signals that 
the designer has chosen to regulate, and that reflect the performance of the system. 
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Figure 3: Nominal controller modified by Q. 


Q is connected to the nominal controller in such a manner (according to the equa- 
tions above) that it sees no feedback. That is, the transfer function from v to e is zero. 
Thus, if Q is stable, it cannot affect the stability of the closed-loop system. By allowing 
Q to range over all stable transfer functions, the combination of the nominal controller 
and Q span the space of all possible stabilizing compensators. 

This ^-parameterization theory can be used in a design procedure by searching for a 
suitable Q p arame ter within a finite subset of the stable transfer functi ons . This method 
has been implemented by Boyd [12] in a computer code called QDES. The Q parameter 
is represented by 

G = X v ‘Qi (21) 

i 

a linear combination of a finite set of fixed, stable maps Thus Q is restricted to a set 
of finite impulse response (FIR) filters. In QDES, numerical optimization is then used to 
search for a set of design variables v that minimize the objective function while satisfying 
the constraints. 

The Q -design approach offers several attractive features. It allows the user to 
directly specify an objective function and strict equality and inequality constraints on a 
wide variety of closed-loop system characteristics. (See Table 2.) In addition, it is 
capable of producing a wide variety of controllers, including those obtainable with LQG 
[15] andH„ [16] methods as subsets. 
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■ Transfer function 2-norm (LQG) 

- Power spectral density of a transfer function at specified frequencies 

> Overshoot and undershoot for a step input 
• Stability margins 

> Infinity-norm of impulse and step responses 

> Transfer function infinity-norm (H-Infinity) 


Limits on impulse and step responses of states and controls 
■ Limits on step response overshoot and undershoot 

• Limits on transfer function infinity-norm 

< Limits on step and impulse response infinity-norms 

• Limits on stability margins 


Table 2: Sample QDES control specifications. 


What makes this approach tractable is the fact that, excluding the structural design 
variables, the resulting optimization problem is a convex one. This is a great benefit; 
even though the number of design variables may be very large, the problem remains rela- 
tively simple and a solution will be found if it exists. 

Admittedly, Q-design has some drawbacks. Because it is based on numerical opti- 
mization, it can require a significant amount of computation. In addition, the resulting 
controllers are of very high order. Nevertheless, the method does show promise. In our 
experience with it so far, the computational effort has been large but not prohibitive, and 
the high order controllers have been amenable to reduction techniques. 

Demonstration Problem 

A representative problem was selected on which to demonstrate the integrated design 
methods and software developed for this task. This example structure, known as the 
Earth Pointing System (EPS) [13] is shown in Figure 4. It is derived from the Ford Earth 
Observation Sciences geostationary platform and consists of a truss-type bus with two 
flexible antennas. The bus is approximately 25 meters in length and the antennas are 15 
and 7.5 meters in diameter. 
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Figure 4: The Baseline EPS Structural Model 


The total mass of the baseline structure is 1027.95 kg, including 150 kg of actuator 
mass and 548.32 kg of nonstructural mass on the antennas. The truss has a 3.0 meter 
cross section and consists of 135 graphite/epoxy tubes. The first flexible mode of the 
structure mostly involves rotation of the antennas and has a frequency of 0.24 Hz. Shapes 
and frequencies of the first 16 flexible normal modes are shown in Appendix A. Many of 
these modes involve local deflections within the antennas with insignificant motion of the 
truss. 

The design problem considered here is one of active vibration suppression. The 
flexible structure has no articulated elements; the controller is required to regulate the 
structure about some nominal attitude while enhancing the stability of the rigid body 
modes. A set of reaction wheel actuators is located near the center of gravity of the base- 
line system and is used to control the pointing errors of the antennas due to the flexible 
response of the structure. Angular rate sensors are colocated with the reaction wheel 
actuators, and linear accelerometers are included at the center of each antenna. 

Dynamic excitation of the structure is considered from a thruster located on the main 
bus module. To expedite the present study, no other loads were considered. The methods 
developed, however, do not restrict the number or characteristics of the external loads. In 
fact, satisfactory practical designs would seem to require many load cases, in both the 
time and the frequency domains. 

A total of 19 structural design variables were chosen for the EPS structure, including 
bar radii for various groups of bars, the width and depth of the truss at several points 
along its length, and the size and orientation of the antenna supports. These design vari- 
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ables are described in Table 3. Figure 5 shows an example of the structure for which 
arbitrary values have been chosen for each design variable. 


Name 

Description 

Baseline 

Value 

C_WIDTH 

Width of the miss at its center. 

■9 

C_DEPTH 

Depth of the truss at its center. 

mSm 

L_WIDTH 

Width of the truss at the large antenna end. 

3.0 m 

L_DEPTH 

Depth of the truss at the large antenna end. 

3.0 m 

SMALL_Y_1 

Value of the ^-coordinate of the truss at the 
end nearest and on the same side as the small 
antenna. 

1.0 m 

SMALL_Y_2 

Value of the y-coordinate of the truss at the 
end nearest and on the opposite side as the 
small antenna. 

1.5 m 

S_DEPTH 

Depth of the truss at the small antenna end. 

0.5 m 

ALPHA 

Angle between the truss axis and the plane of 
the antennas. 

0.0 

LONGERON_R 

Radius of the longerons. 

0.0255 m 

BATTEN_R 

Radius of the battens. 

0.0255 m 

DIAG_EXT_R 

Radius of the external diagonals. 

0.0255 m 

DIAG_INT_R 

Radius of the internal diagonals. 

0.0255 m 

L_SUP_R 

Radius of the large antenna supports. 

0.0255 m 

S_SUP_R 

Radius of the small antenna supports. 

0.0255 m 

ACT_SUP_R 

Radius of the actuator supports. 

0.0255 m 

ACT_MASS 

Actuator mass. 

150.0 kg 


Table 3: Design variables for the EPS structure. 
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L — — — 

Figure 5: Structure with arbitrary changes in the design variables. 

The EPS equations of motion that were used for the control design contained the first 
16 flexible modes, ranging from 0.24 Hz to 3.53 Hz. A sample frequency of 20 Hz was 
used, with 20 Q taps and 500 samples. The sample rate is sufficient to capture at least 5 
points in every cycle of the highest frequency mode, and the 500 samples used for the 
open-loop impulse responses is sufficient to capture 6 full cycles of the lowest frequency. 

Figure 6 illustrates the locations in the baseline model relevant to the control design. 
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Rotations about the x, y and z axes are designated as 6, <j)> and \ff, respectively. The 
actuators are reaction wheels located at node #42. The torques exerted by the actuators 
are designated as/ 042 ,/^ 2 , and/^42- A thruster, which represents a disturbance for the 
control task, is located at node #45 and exerts a thrust / Z 45> tn the positive z -direction. 
The measured outputs are the three components of angular rate at node #42 and the linear 
accelerations in the x and z -directions at nodes #26 and #71. The control objective is to 
regulate the antenna angles Y^26> ^ 71 , and y/j\ using the reaction wheel actuators. 

The initial system configuration provided by the structural model consists of 1 
exogenous input, 3 controller inputs, 4 regulated outputs, and 7 measured outputs. This 
configuration is depicted in Figure 7. 
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To accommodate the design of a robust system, several additional inputs and outputs 
were added to the system. As shown in Figure 8, actuator noise and sensor noise were 
added. Also, the exogenous output was augmented with the noisy actuator signal and the 
angular rates measured at the reaction wheels. This configuration now allows for a wide 
variety of control specifications, including robustness and noise sensitivity. 



Findings 

The COSTAR implementation of the integrated design method described previously 
was applied to the EPS design problem. This was done in two ways: (1) with all the 
structural design variables held fixed, and (2) with all but one of the structural design 
variables held fixed. The first analysis amounts to an optimal control design for the 
baseline structure, which is useful for evaluating the (2-design methodology. The second 
analysis is an integrated design optimization of the controlled EPS structure. Because the 
analytical sensitivity computations were incompletely implemented in COSTAR at the 
time, the latter analysis considered only a single concurrent structural design variable. 
Although this is a severe limitation in the optimization problem, the integrated design 
capability of COSTAR is nevertheless demonstrated adequately. 

Optimal Control Using ^-Design 

By fixing the structural design variables, COSTAR is essentially identical to the 
QDES control design code. The control design variables shown in Equation (21) are 
adjusted to minimize the objective function while satisfying the constraints. This control 
design is for a specific disturbance — a 100 N-s impulsive force in the z-direction by the 

thruster at node #45. 
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The objective function chosen for this demonstration contained control effort, point- 
ing error, and stability margin components. Specifically, this objective was calculated as 
the sum of 1(P times the mean square torque, 10 16 times the mean square pointing error, 
and 10 4 times the inverse of the minimum Nyquist distancef . (All quantities are mea- 
sured in SI units.) These costs were chosen for a disturbance environment where the ratio 
of thrust disturbance to sensor noise is approximately 10,000. 

Constraints were applied to the antenna pointing error, the actuator torque levels, and 
robustness to actuator noise. The maximum allowable pointing error of either antenna 
was 0.01°. The torque applied by the reaction wheel actuators was limited to at most 
1000 N m. The minimum Nyquist distance was required to be at least 0.5 S" 1 . 

These objective and constraint functions, along with the impulse responses of the 
baseline system, were given as input to QDES. The result was an optimal controller that 
minimized the costs without violating the imposed constraints. As mentioned earlier, 
QDES often results in high order controllers, and this design was no exception. For the 
36 state baseline system, QDES found a 146 state controller. It was demonstrated using 
the controller order reduction technique of frequency weighted internal balancing [14], 
that the controller could be reduced to 30th order with minimal impact upon performance. 

To the closed-loop system thus obtained, a 1 N s impulsive disturbance in the 
z-direction was applied by the thruster at node #45. This is precisely the disturbance for 
which the controller was designed, except that its magnitude is scaled by 100. The result- 
ing response and the response of the open-loop system are depicted in Figures 9-18. 

Figure 9 depicts the pointing error about the x -axis at the small antenna, £71, for this 
impulsive disturbance. Recall that the control specifications included a constraint on 
pointing error of 0.01° for a 100 N s impulsive input. This translates into a 1.7xl0~ 6 
radian limit for a unit impulse, which is shown in the figure. One can see that the open- 
loop response clearly violates this constraint. It has a peak response of —6.5x10"^ radians. 
On the other hand, the closed-loop response goes right up to the constraint in about 0.5 
seconds, but never violates the constraint. Of the four pointing angles, fyi is the most 
sensitive to the thruster disturbance. The other pointing angle responses, shown in 
Figures 10-12, do not violate the pointing constraints even for the open-loop case. The 
response of y/26, the pointing error about the z -axis depicted in Figure 12, is roughly 
three orders of magnitude less than that of #71. 

Figures 13-15 illustrate the effect of the impulse on the angular rates measured at the 
actuators. The closed-loop response exhibits greatly improved dynamic behavior. No 
constraints were imposed upon these angular rates. 
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Finally, Figures 16-18 depict the response of the actuators to the impulsive distur- 
bance. The actuators were constrained to have a 10 N m peak torque level. For the base- 
line structure, however, none of the actuators approached this limit. A peak torque level 
of approximately 2 N m was seen in actuators and/^2- During the combined opti- 
mization, however, it is entirely possible that the structure might evolve into one where 
this peak torque limit could be reached. Placing this limit prevents the optimization pro- 
cedure from designing a system that would require excessive actuator authority. 
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Figure 9: Response of (h\ to a unit/ Z 45 impulse. 


t The minimum Nyquist distance is the minimum distance between the Nyquist curve 
and the “-1 point” (the point of instability) in the complex plane. 
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Figure 10: Response of yj\ to a unit / Z 45 impulse. 



Figure 11: Response of 026 to a unit/^s impulse. 
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Figure 12 : Response of y/fo to a unites impulse. 



Figure 13 : Response of 64 2 rate to a unit .£45 impulse. 
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Figure 15: Response of y /42 to a unit/ Z 4s impulse. 
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Figure 18: Response of fy 42 to a unites impulse. 

Combined Optimization 

The results presented above pertain to a COSTAR optimization in which only control 
design variables were considered. By extending that analysis to include structural design 
variables, the integrated design capability of COSTAR can be demonstrated. As men- 
tioned previously, the analytical sensitivity computations were incompletely implemented 
in COSTAR, forcing us to consider only one structural design variable at a time. We per- 
formed two such analyses, including the width and depth of the center section of the EPS 
truss. (These are the variables CJWIDTH and C_DEPTH described in Table 3.) 

For this integrated design problem, the objective function described in the previ- 
ous section was augmented to include a term proportional to the total mass of the struc- 
ture. In many cases, this mass term may adequately account for the cost of the materials 
and the cost of transporting the structure to orbit. It cannot account, however, for many 
other costs such as manufacturing and assembling the structure. As with any structural 
optimization procedure, the results are only as good as the objective function. Although 
this demonstration used a very simple objective function, the only objective function limi- 
tation in COSTAR is that it must be twice-continuously differentiable. 

In addition to the control constraints described previously, simple side constraints 
were applied to the structural design variables. These constraints specified commonsense 
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restrictions such as that the truss cross section must have a reasonable width and depth. 
Far more complicated constraints (on stresses, strains, etc.) could be included if desired. 

In Figure 19, the objective value is shown as a function of the center section truss 
depth. Each point on the curve represents a particular design analyzed by COSTAR in 
which an optimal control design is obtained. The baseline EPS structure has a depth of 
3.0 m; the figure shows that its optimum value (for the objective and constraint functions 
specified and for the particular finite-dimensional space from which Q is chosen) is 
approximately 2.6 m. Significant contributions to the objective function value come from 
the total mass and mean square pointing error terms. Because the mass varies by only 7% 
over the range of truss widths shown, it does not play a significant role in determining the 
location of the optimal design. Nevertheless, incorporation of the structural design vari- 
able clearly has an important effect on the closed-loop behavior and on the performance 
of the resulting design. 



Depth (meters) 

Figure 19: Objective vs. truss depth. 

Figure 20 shows the combined control/structure optimum configuration obtained. At 
2.6 m, the center section depth has decreased from the baseline EPS design — a change 
that significantly decreases the bending stiffness of the truss. This narrowing of the cen- 
ter section obviously decreases the structural cost (total mass) slightly, but it also 
decreases the combined costs associated with control energy and performance. 
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Figure 20: Optimized structure considering control and truss width. 


A similar analysis was performed with the center section width, rather than depth, as 
the structural design variable. A plot of the objective value verses this width is shown in 
Figure 21. Here again, each point represents a system for which an optimal controller has 
been designed with COSTAR. The baseline EPS structure has a depth of 3.0 m; the fig- 
ure shows that its optimum value is approximately 1 .4 m. 



Figure 21: Objective vs. truss width. 

Because the optimal control problem is convex, we know that each point in Figure 21 
represents a globally optimal design for that particular structure. As the figure demon- 
strates, however, the optimization problem is no longer convex when structural design 
variables are introduced. Several local minima are seen to exist. Clearly, finding a global 
optimum is not generally assured for the combined problem. 
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Figure 22 shows the combined control/structure optimum configuration considering 
the truss width structural variable. Note that, as with the center section depth, the opti- 
mum center section width is also less than that of the baseline EPS structure — again a 
change that significantly decreases the bending stiffness of the structure. 



Figure 22: Optimized Structure Considering Control and Structural Depth 

Conclusions 

In conducting a simultaneous optimization of both structural and control parameters, 
it is possible to investigate design configurations which a separate control and structure 
optimization would miss. Thus, the overall approach to combined control/structure opti- 
mization presented in this report holds promise for generating non-intuitive optimal 
designs. In particular, the optimization of the EPS structure described in this report found 
that it was optimal to decrease the size of the center truss section. This represents a 
reduction in structural stiffness - a result which might run counter to intuition. 

In the process of developing this combined control/structure optimization, a method 
for determining structural sensitivities was developed. This method of analytically calcu- 
lating structural sensitivities represents a significant improvement over other methods 
involving either finite differencing or the differentiation of eigenvectors. 

The control aspect of the optimization was accomplished using the relatively new 
control synthesis tool called Q-design. As implemented by the program QDES, this 
design method has been demonstrated to work well in the combined optimization context. 
It has proven to be a very flexible tool for making control design specifications and syn- 
thesizing controllers that meet these specifications. 
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With this increased flexibility also comes increased computational cost. In addition, 
controllers developed using Q-design tend to be of very high order — much higher than 
that of the systems they are controlling. On the other hand, it has been demonstrated that 
the controller order reduction technique of frequency weighted balancing can greatly 
reduce controller order. Moreover, the use of Q-design in the combined optimization 
does not preclude the use of other control design tools for the optimized structure; the 
benefits of Q-design would have then already been exploited by designing a system for 
which it is theoretically possible to achieve the desired performance. 

While optimization over the control design variables is a convex problem, it is not 
convex over the structural parameters. Consequently, the combined control/structure 
optimization will most likely have local minima, and the process of attempting to find a 
design that achieves a global minimum will require multiple starting points. 


Areas for Future Research 

At this point, the COSTAR optimization procedure requires a great deal of supervi- 
sion and designer input. This is related primarily to facilitating communication between 
the various software components. One of the primary goals of future work would be to 
streamline the optimization procedure. Ideally, once the structural model and its design 
variables have been defined, and once the objective and constraint functions have been 
specified, the optimizer should be able to run "hands free until an optimal design is 
found. This streamlining of COSTAR would also facilitate the inclusion of a larger num- 
ber of structural design variables. Present designs have only been optimized over one 
structural variable and all control variables. 

Some of the subprocesses in the optimization procedure have room for much 
improvement. In particular, the implementation of analytic design sensitivity computa- 
tions should be completed in order to facilitate the inclusion of a larger number of struc- 
tural design variables. Another area for research is the eigenvalue extraction method, 
which might make better use of the eigenvalues calculated for previous structures, since 
from one design iteration to the next the structure does not change a great deal. The 
effects of modal truncation, as they relate to this simultaneous optimization problem, 
should be studied further. 

Finally, it would be advantageous to develop a means of reducing the computation 
time associated with the control part of the optimization. Some research has been done in 
this area by Kosut and Kabuli, and is included as Appendix B in this report. This 
approach involves a pre-filtering of the (7-filter in an attempt to reduce the number of taps 
necessary to obtain a good control design. Another approach would be to incorporate a 
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nominal controller for the baseline system. This would reduce the number of samples 
required to characterize the impulse responses, and would also reduce the number of taps 
required in the Q-filttr. 


References 

[1] Vos, R.G.; Beste, D.L.; and Gregg, J., “Integrated Analysis Capability (IAC) User 
Manual (Level 2.5)”, Boeing Aerospace & Electronics, Seattle, Washington, 
1989. 

[2] Messac, A., “Optimal Simultaneous Structural and Control Design of Large Space 
Structures”, Ph.D. thesis, Massachusetts Institute of Technology, Cambridge, 
Massachusetts, November 1985. 

[3] Khot, N.S.; Venkayya, V.B.; and Eastep, F.E., “Optimal Structural Modifications 
to Enhance the Active Vibration Control of Flexible Structures”, AIAA Journal, 
Vol. 24, Apgust 1986, pp. 1368-74. 

[4] Miller, D.F. and Shim, J., “Gradient Based Combined Structural and Control 
Optimization”, Journal of Guidance, Control, and Dynamics, Vol. 10, May-June 
1987, pp. 291-8. 

[5] Onoda, J. and Haftka, R.T., “An Approach to Structure/Control Simultaneous 
Optimization for Large Flexible Spacecraft”, AIAA Journal , Vol. 25, August 
1987, pp. 1133-8. 

[6] Morrison, S.K.; Ye, Y.; Gregory, C.F.,Jr.; Kosut, R.L.; and Regelbrugge, M.E., 
“Integrated Structural/Controller Optimization for Large Space Structures”, AIAA 
Guidance, Navigation, and Control Conference , Minneapolis, Minnesota, 1988. 

[7] Milman, M.; Salama, M.; Scheid, R.; Bruno, R.; and Gibson, J., “Integrated 

Control-Structure Design: A Multiobjective Approach”, JPL D-6767, Jet 

Propulsion Laboratory, California Institute of Technology, January, 1990. 

[8] PATRAN Users’ Guide, PDA Engineering, Santa Ana, California, September 
1986. 

[9] MACSYMA Reference Manual , Document No. SMI0500030.013, Symbolics, Inc., 
Burlington, Massachusetts, November 1988. 


32 



[10] Dailey, R.L., “Eigenvector Derivatives with Repeated Eigenvalues”, AIAA 
Journal, Vol. 27, No. 4, April 1989, pp. 486-91. 

[11] Youla, D.C.; Jabr, H.A.; and Bongiomo, J.J., “Modern Wiener-Hopf Design of 
Optimal Controllers — Part II: The Multivariable Case,” IEEE Transactions on 
Automatic Control, Vol. 21, pp. 319-38, June 1976. 

[12] Boyd, S.P.; Balakrishnan, V.; Barratt, C.H.; Khraishi, N.M.; Li, X.M.; Meyer, 
D.G.; and Norman, S.A., “A New CAD Method and Associated Architectures for 
Linear Controllers,” Technical Report No. L-104-86-1, Stanford University, 
December 1986. 

[13] Belvin, K., “EPS Structure Description”, NASA Langley Research Center, June 1, 
1989. 

[14] Gregory, C.Z., “Reduction of Large Flexible Spacecraft Models Using Internal 
Balancing Theory”, AIAA Guidance and Control Conference, Gatlinberg, 
Tennessee, August 1983. 

[15] Skelton, R.E., Dynamic Systems Control: Linear Systems Analysis and Synthesis, 
John Wiley and Sons, 1988. 

[16] Francis, B.A. and Doyle, J.C., “Linear Control Theory with an //« Optimality 
Criterion”, SIAM J. Cont. Optim., Vol. 25, No. 4, July 1987, pp. 815^14. 


33 


Appendix A 

This appendix contains plots showing the mode shapes of the first 16 flexible modes 
of the baseline EPS structure. These are screen copies from an interactive PATRAN 
session. 
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1 Introduction 


Consider a linear time-invariant finite-dimensional discrete-time plant described by the 
transfer function P (see Figure 1); w , z , u and y denote the exogenous inputs, regu- 
lated outputs, actuator inputs and sensor outputs, respectively. Typically, P represents 
the augmented plant', the original plant Pyu (the block of P representing the transfer 
function from u to y) is augmented with stable blocks used as weights when taking into 
account the exogenous inputs and regulated outputs. In general, all that is required is that 
the unstable poles of P are in P yu . With this assumption (which is in fact a necessary 
and sufficient condition for internal stabilization of the closed-loop system in Figure 1 by 
dynamic compensation) the set of all compensators C that stabilize P in Figure 1 is 
the set of all compensators that stabilize P yu . 


w 

u 



z 

y 


Figure 1: Closed-loop system 

For a given lightly-damped flexible space structure with position and/or rate measure- 
ments, the transfer function P yu is stable (poles in the open unit disk) and strictly-proper 
( P yu ( oo) = 0 ). Hence, the set of all proper compensators C such that the closed-loop 
system in Figure 1 is stable is given by 

{ Q(I — PyuQ) -1 | Q is proper and stable } . (1) 

Using the parametrization in (1), the set of all achievable stable closed-loop transfer func- 
tions H zw is given by 

{ Pzw ~ PzuQPym \ Q is proper and stable } . (2) 
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2 Problem Description 


Once the complete solution to all achievable closed-loop transfer functions is parametrized, 
one can conceivably improve (if possible) the closed-loop performance by reducing an 
objective function of H zw over the parameter Q (from now on referred to as the Q- design). 
However the set of all realizable stable transfer functions can not be parametrized by a finite 
number of real parameters. Nevertheless, by increasing the number of real parameters, a 
sequence of parameter optimizations can be performed until the desired improvement in 
performance is achieved (as in QDES). Note that for an rc; input n 0 output plant P r , , 
an .Nth order FIR approximation requires N • rii ■ n 0 parameters. 

The problem is to find a “suitable” class of stable transfer functions {Q g ), e for the 
specific plant at hand such that N is “small” and minimization over q € IR N achieves 
the desired performance improvement. 

Provided that a “suitable” parametrization of stable Q is available, the strong points 
of Q-design can be summarized as follows: 

• As long as the design problem is approximated as a parameter optimization prob- 
lem, the objective function can include both time- and frequency-domain constraints 
(choosing the objective function convex in the parameters definitely helps). Such a 
mix is not suitable in the conventional PC2~ a-nd 7Yoo-design methods. 

• If the input-output description is already given by impulse responses, there is no 
need for a model identification until the compensator is to be implemented. Since 
the feedback stabilization problem is already taken care of, each objective function 
evaluation need not be done after a stabilizing feedback law is calculated using con- 
ventional design methods, relying on the identified model. 

• All three of the transfer functions in (2) (namely, P zw , P zu and P yw ) depend on the 
open-loop plant (the weights introduced during augmenting P yu are stable cascaded 
blocks which do not introduce any feedback paths). Hence if the plant is also to be 
parameter optimized, the plant parameters are simply appended to those of Q . 

In the rest of this report, we focus on a case study and a simplified SISO model of 
the EPS (Earth Pointing Satellite). As a benchmark comparison, Q-design is compared 
against ^-design. 


52 


3 Method used for SISO P yu 

Suppose that there is one actuator and one sensor. Let the transfer matrices Ti and T 2 
be defined as follows: 

T\ := P zw 

T2 '.= PzuPyw • 

The following steps are taken to compute the optimal (in the ?f 2 -norm sense, over a 
parametrization Q q : q £ 3R N ) norm of the transfer matrix 

H zw — T\ — Q q T2 : 


1. Determine the number of samples N sa . mples ; get the impulse response sequences for 
the transfer matrices T\ and T 2 . 

2. Choose a numerator polynomial n and a strictly-Hurwitz denominator polynomial 
d of the same order. Determine the number of taps IV taps . Set 


N, 


Q« = ? Y. 1 ; 

° «=i 

we choose the FIR term strictly-proper since the optimal Q for a proper plant is 
strictly-proper. 

3. Filter the pulse sequences of T 2 with ^ : 


taps 


n 


T2- . 


4. Determine a solution q of 

^ taps 

argmin||Ti - T 2 ^ z ~' IU 

from the minimum norm least squares solution of 

9i 

y = A \ 

9iV, 


taps 

where y and A can be obtained by reordering the entries of Ti and T 2 . To 
find the solution, there is no need construct y and A completely; in our cases 
w , z £ IR 2 and iV samp i es = 5000 and 10000 . Instead, form A T A £ iR-^tapsxNtaps 
and A T y £ IR‘^ taps . Get the minimum norm least-squares solution 

q - (A T A)\(A T y) . 

Compute the least-squares error ||y — Aq\\ 2 without constructing A . 
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5. Plug the compensator 


C = Q#(l - PyuQi)-' 


in the feedback loop. 

6. Compute the ?f 2 -norm analytically from the state-space description of the closed- 
loop map H zw and check against step 4 . 
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4 Case Study 


Let G denote the zero-order hold equivalent (sampled at 50 Hz 1 ) of the transfer function 
3 i + ois+i ■ Let P in Figure 1 be given by 


’ Zl ' 


100G 0 100G ' 


' Wi " 


=z 

0 0 0.1 


w 2 

. y . 


G 0.01 G 


u 


(3) 


Discrete-time ^-optimal design gives | \H ZW \\ 2 , op t = 0.8225 . 
We chose N samp i es = 5000 ( 5 time constants). Let 


n = d = 1 . 


For a given iV taps , let 


Q< = 


-^taps 

• 


i=l 


For different number of taps, H-ff^lkEXE (the computed least-squares error ||y — Aq\\ 7 
using the executables) and 1 1 ik™ I ksYS ( W 2 - norm from the system matrix description of 
the closed-loop map H zw obtained by plugging in the compensator Q q { 1 — P yu Q q )~ l •) 
are listed in Table 1 . 


^taps 

ll-ffi^lb.EXE 

ll-fl^tulksYS 

l 

31.5830 

31.5839 

5 

2.1230 

2.1241 

10 

1.6210 

1.6215 

20 

0.8513 

0.8514 


Table 1: Optimal FIR approximations for different number of taps;optimal 7Y 2 -norm is 
0.8225. 

Note that column 2 in Table 1 is computed from the 5000 samples of the impulse 
responses of P zw , P zu and P yw . Column 3 is obtained from the closed-loop system 
in Figure 1, where the plant P is represented by the state- space description. The close 
results are due to considerably high sampling rate and long pulse sequence used in the 
computations. For a larger scale problem, one should be cautious in assigning the sampling 
rate and sample size, trading off computational burden versus aliasing. 

1 Continuous-time 7f 2 -optimal design for the plant description where G in (3) is replaced by jy^YI+I 
results in a compensator eigenvalue magnitude of 41.56 rad/s . The sampling rate was chosen approximately 
7.5 times faster to get the discrete-time optimal design performance similar to the continuous-time frequency 
plots. 
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In order to reduce the number of parameters ( N Ups ), an initial check was made by 
assigning 

n — z 

d = n ( z -pd 

»=i 

j 2—/ 9* z ’ 

a .=i 

where the preassigned p; ’s denote the optimal Q pole locations. Solving for this 4 
parameter Q representation ( Q q ), we obtained 

||H 2ti ,|kEXE = 0-8306 
||^|| 2 ,sys = 0.8306 . 

Clearly, this 4 parameter representation of Q achieves a better result than the 20 pa- 
rameter FIR representation in Table 1 . However, one does not have access to the optimal 
pole locations to start with. 

A couple of other d assignments were done as follows: for a fixed radius r , the four 
poles were spread as re ±j9k where 8 k £ [0 , 0 max ] . For \6i 6 2 ] = [0.1 0.5] , we obtained 
||#z«;|| 2 ,exe = ||#zu,|ksYS = 1 0645, 0.8969 for r = 0.95 and 0.8, respectively. Both of 
these trials did better than a straightforward 4th order FIR approximation (2.89) . 

The motivation is to preassign a structure for Q as in Section 3 (typically a stable 
factor of the inverse of the plant) and fine tune with the FIR factor, hopefully over a 
smaller number of parameters. If the openloop plant is to be designed taking into account 
the performance objective desired from a feedback law, this preassigned structure will be 
at most in the order of the plant and depend only on the plant parameters. Note that 
an identification approach in assigning the numerator and denominator coefficients of 0 
may not be useful due to the constraint that the parameter Q is stable. 
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5 EPS Model 


In order to try out the idea involving the plant inversion, an SISO approximation of the 
original 32-state EPS model is obtained. The procedure is explained below. 

Consider the 32-state state space description of EPS: 


X — Ax ”1” -Bthrdthr T Hact^act 


yiosi 


^los^ 


yios2 


^thr 

:= /z45 


^act 

:= [fe 42 

f<t> 42 /tM2J 

yiosi 

:= [$26 

fcef 

yios? 

'•= \$71 

M T ■ 


T 


(4) 


There are 16 modes ranging from 0.2423 Hz to 3.5317 Hz . The damping ratio is 

£ = 0.02 . 

The four singular-value plots of the openloop plant in (4) are shown in Figure 2 . 


5.1 Simplified EPS Model 

We used the following 3 mode reduced model in (5) for the EPS : 

X = 

yios — 
u :~ 

yios • “ 

obtained by the following MATRIXx commands: 

SORIG32=[A BTHR BACT;CLOS 0*EYE(4)]; 

[SBAL SIG T]=BALANCE(SORIG32([1:10 35], [1:10 34]),10); 

[AA BB CC DD]=SPLIT(SBAL,10); 

BTHR=T\SORIG32(1:10,33); 

ABAR = AA(1:6,1:6); 

BBARTHR - BTHR(1:6) ; 

BBARACT = BB(1:6); 

CBARLOS = CC(:,1:6); 

The discrete- time model used for the EPS model is the zero-order hold equivalent 
(sampled at 300 Hz) of the 6-state plant description in (5). 


Ax T Rthr^thr T 
fe 42 

$71 


(5) 
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The magnitude plots for the discrete-time 6-state EPS model are shown in Figure 3. 

Compared with plots C and D in Figure 2 , plots A and B in Figure 3 match the 
frequency responses up to approximately 6 rad/s . 
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6 Closed-loop Performance of SISO 6-State Discrete- 
Time EPS Model 

Let G denote the zero-order hold equivalent (sampled at 300 Hz) of the 6-state plant 
description in (5) where 

yios = Cndt.hr + G 2 u • 

Consider the augmented plant model 

' 100G1 0 ioog 2 i r wi ' 

0 0 0.1 w 2 . (6) 

Gi -0.01 G 2 J [ u 

Apply the feedback law 

u = —Q q { 1 — G 2 Q q ) 1 y 

as in Figure 1 , where Q q is to be determined. 

As a reference performance measure, discrete-time ^-optimal compensator is de- 
signed; the closed-loop yields ||H’ ztt ,||2 |0 pt = 2.284 . The magnitude plots for the four 
entries of the Tib-optimal H zw axe shown in Figure 4. Comparing D in Figure 4 and B in 
Figure 3, the optimal Q magnitude plot is the negative of the P yu plot (over the 0.3 rad 
range in Figure 3; taking into account the 60 dB shift due to the weights 0.1 and 0.01). 

Sample size is chosen as N samp \ es = 10, 000 (approximately 1 time constant at 300 Hz). 
As in the case study in Section 4 , FIR approximations are made for different number of 

taps ( Q q = Qi z ~' )• Results are listed in Table 2: 



Table 2: Optimal FIR approximations for different number of taps;optimal Tf 2 -norm is 
2.2840. 

Note that, unlike Table 1, H-ff^lh.EXE and H-ff^l^sYS (see section 4) are not close, 
due to 1 time constant truncation of the original pulse sequence. Despite large N taps there 
is not a considerable improvement in the closed-loop 7i 2 - norm. 

Instead, we chose 

Q q = 

d i = 1 

where n has all six poles of G 2 and d has all four stable zeros of G 2 and two at zero 
(the two zeros at zero are chosen to make 4 relative degree zero; no specific reason for 
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assigning them at zero). Hence ^ represents the stable factor of the inverse of Gi and 
the four FIR parameters are introduced to fine tune. Calculations show that 


| \H zw | |2,exe = 7.7974 

||fl*w|ksYS = 3.2685 (optimal is 2.284) . 

The magnitude plots of the four entries of the final run H zw with Q q 
in Figure 5. 


as in (7) axe shown 


62 


c .0001 


.001 .01 .1 
rad 


! D .0001 .001 .01 .1 1 

rad 

Figure 4: 


The magnitude plots for the four entries of the 7f 2 -optimal H zw for the plant description 
in (6): A - w x z r , B - w 2 ■-» z x , C - wi z 2 , D - w 2 »-+ z 2 . 






The magnitude plots of the four entries of H zw for the plant description in (6) and Q 
in (7): A - W\ t— ► z\ , B - w-i t— * z\ , C - u>i *— ► z 2 , D - wi •— »• zz . 
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